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We investigate contributions to the extragalactic gamma-ray background (EGB) due to neutralino 
dark matter (DM) pair-annihilation into photons, from DM density enhancements (minispikes) sur- 
rounding intermediate-mass black holes (IMBHs). We focus on two IMBH formation scenarios; our 
conservative scenario where IMBHs are remnants of Population-Ill stars, and our optimistic scenario 
here IMBHs are formed in protogalactic disks. In both scenarios, their formation in pregalactic halos 
at high redshift lead to the formation of minispikes that are bright sources of gamma-ray photons. 
Taking into account minispike depletion processes, we only sum contributions from a cosmological 
distribution of IMBHs with maintained minispikes. Our conservative scenario (BH mass 1O 2 M 
with a r -3 / 2 minispike) predicts gamma-ray fluxes that are an order larger than the equivalent flux, 
using the same DM parameters (mass 100 GeV and annihilation cross-section 3 x 10~ 26 cm 3 s _1 ), 
from the host halo without IMBH minispikes. Our optimistic scenario (BH mass 1O 5 M with a 
r~ 7//3 minispike) predicts fluxes that are three orders larger, that can reach current EGB observa- 
tions taken by EGRET (DM parameters as above). This fact may serve interesting consequences for 
constraining DM parameters and elucidating the true nature of IMBHs. Additionally, we determine 
the spectra of DM annihilation into monochromatic gamma-rays, and show that its flux can be 
within observational range of GLAST, providing a potential 'smoking-gun' signature of DM. 



I. INTRODUCTION 

Despite compelling indirect evidence, from galactic to 
cosmological scales, the fundamental nature of the dom- 
inant non-baryonic component in the matter density of 
the universe (dark matter, hereafter DM) remains un- 
known. Intriguingly, extended models of particle physics 
independently provide us with a host of particle candi- 
dates for this as yet unknown matter, of which the most 
popular is the supersymmetric neutralino (see reviews 
PIE LJ for details). Upgrades of underground direct de- 
tectors looking for scattering of DM particles from nuclei, 
together with future neutrino, antimatter, and gamma- 
ray detectors looking for products of DM annihilation, 
will dramatically enhance our chances of understanding 
the true nature of DM. In particular, the forthcoming 
launch of the Gamma Ray Large Area Space Telescope 
(GLAST) || and numerous ground based Atmospheric 
Cerenkov Telescopes make indirect gamma-ray search es- 
pecially promising. 

Since the DM annihilation rate scales as the DM den- 
sity squared, there is great advantage in observing areas 
where the DM density is believed to be high. The galactic 
centre (GC) is the immediate choice, and indeed strong 
gamma-ray emission has been observed and its nature 
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and origin have been investigated by many researchers 
d B HE H Elm IH El However, the DM density in 
the GC is highly uncertain, making accurate predictions 
difficult. For example, DM enhancements called 'spikes' 
can form during the formation of a central supermassive- 
BH (SMBH) [3, but it can also be depleted by various 
processes by varying degrees m> El O E3- I n addi- 
tion, nearby astrophysical (non-DM) gamma-ray sources 
make a potential DM detection imp ossible for all but a 
narrow range of DM parameters |l9l . Intermediate-mass 
black holes (IMBHs, see e.g. |20j) provide an alterna- 
tive source that may work positively for DM detection. 
Bertone et al. [2l| recently investigated the possibility of 
detecting a 'smoking gun' gamma-ray signature of DM 
using IMBHs in the Milky- Way as point sources. They 
showed that IMBH formation increases the DM density 
in its vicinity to produce a 'minispike', and also that DM 
enhancement depletion processes are generally less signif- 
icant for IMBHs due to their roughly spherical distribu- 
tion about the GC. They conclude that under optimistic 
circumstances, the Energetic Gamma Ray Experimental 
Telescope (EGRET) may have already seen a few of the 
IMBH minispikes as unidentified sources. 

Another avenue of indirect DM search is via the extra- 
galactic gamma-ray background (hereafter EGB) mea- 
sured over a wide energy range j^, . The origin of 
this background is currently unknown, and it has been 
speculated that DM annihilation gamma-rays from cos- 
mological distributions of DM contribute to some degree 
0mmmillllim& Since the DM annihi- 
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lation cross-section is so small, a consideration of DM 
enhancements is crucial for meaningful gamma-ray flux 
predictions. A popular DM enhancement is those at the 
centres of galactic DM halos. However, Ando [2^| has 
recently shown that they are strongly constrained by ob- 
servations of our galaxy. The author assumes universality 
of galactic DM halo profiles, and shows that DM annihi- 
lation cannot significantly contribute to the EGB with- 
out exceeding gamma-ray observations from our GC [2^ . 
The author also points out that this constraint could be 
loosened when one takes DM substructures into account. 

In this paper we argue IMBHs minispikes as a substruc- 
ture in the DM halo, that can lead to enhancements that 
do not conflict with current observations of our galaxy. 
We determine contributions to the EGB from IMBH min- 
ispikes by summing gamma-ray fluxes from all redshifts. 
IMBH minispikes are not expected to greatly suffer from 
depletion processes, but we do take into account BH-BH 
mergers, which are known to strongly deplete minispikes 
and do occur in IMBHs. We also consider a conserva- 
tive case (10 2 Af Q BHs of Population-Ill origin [3^ with 
a r~ 3 / 2 minispike) and an optimistic case (~ 10 5 M Q BHs 
formed in the centres of protogalactic disks [3i| with a 
J--7/ 3 minispike). Our result is that contributions to the 
EGB are increased by 1 — 3 orders in magnitude. In 
particular, our optimistic case predicts fluxes that can 
reach current EGB observations. As this has interesting 
implications for constraining DM parameters and IMBH 
scenarios, we critically assess uncertainties in our calcu- 
lation. We then determine the flux of DM annihilation 
into line gamma-rays, and show that under optimistic 
conditions, it is observable by GLAST. This provides a 
potential 'smoking-gun' signature of DM. 

This paper is structured as follows. In Sec. [H] we in- 
troduce IMBHs, starting with their existence, followed by 
their formation scenarios, and finishing off with a sum- 
mary of recent numerical studies. Then in Sec. II I II we 
develop our calculation frameworks, first for the EGB, 
followed by DM annihilation, then minispike formation, 
moving finally on to our IMBH number density fitting. 
Calculations and results are in Sec. IIVI and discussions 
and conclusions in Sec. In all our calculations we 
adopt the standard flat cosmological constant plus cold 
DM (ACDM) cosmology, with tt M = 0.3, tt A = 0.7, 
h = 0.7, and a s = 0.9. 



II. IMBH: EVIDENCE AND PROPERTIES 

A. Evidence for IMBH 

Clues for the existence of a class of BHs with masses 
heavier than stellar BHs but lighter than supermassive 
black holes (SMBH) have accumulated in recent years. 
We call these BHs intermediate-mass BHs (IMBHs) and 
loosely define their mass range as 20 < M^/Mq < 10 6 . 
We briefly discuss the observational evidences and theo- 
retical motivations for their existence. 



Observationally, studies of objects known as ultra- 
luminous x-ray sources (ULXs, |3J,|33) reveal that they 
may harbour IMBHs. Although most x-ray sources can 
be understood as accretion by compact objects such as 
BHs, there is an upper limit on the luminosity for a given 
BH mass, known as the Eddington limit. For 20Af Q , 
commonly accepted as the upper limit of stellar BHs, 
this limit is ~ 2.8 x 10 39 ergss _1 . ULXs are observed 
to exceed this limit. This phenomena can be explained 
by several mechanisms, including a short-term super- 
Eddington phase |36j| . beaming |37|. or normal accretion 
by an IMBH. Although the debate has not been settled, 
evidence favouring the IMBH mechanism over the other 
two have accumulated in recent years; these come in var- 
ious forms, including spectral analysis pol l38|. evidence 
for a low temperature (~ 0.1 keV) black-body compo- 
nent [3^ , anal ysis of break frequencies of the power den- 
sity spectrum [4(j , and observation of broad Fe lines and 
quasi-periodic oscillations (QPO) [4l(. It seems that at 
least a fraction of the ULXs, in particular the most lu- 
minous ones, are IMBHs. The BH mass inferred from 
ULX observations is of the order ~ 10 3 M Q (assuming 
no beaming and accretion efficiency ~ 10%). Also, from 
their positions in the host galaxies we can deduce an up- 
per limit of < 10 6 M Q , in order not to sink to the galactic 
centre by dynamical friction within a Hubble time |20| . 

Theoretically, the existence of a population of IMBHs 
helps explain the origin of SMBHs. There is as yet no 
definitive SMBH formation scenario, but the discovery 
of quasars at redshifts z <~ 6 in the Sloan Digital Sky 
Survey |H El El |H Ejj suggests that they were al- 
ready formed at high redshifts. Such an early formation 
lends itself to scenarios with massive seed-BHs, frequent 
mergers, and rapid accretion. A hierarchical formation 
scenario starting from massive seed-BHs also helps ex- 
plain the tight correlations observed between the SMBH 
mass and properties of the host galaxy and halo [3^ . The 
natural outcome of this hierarchical scenario is the exis- 
tence of 'wandering' BHs in galactic halos, resulting from 
seed-BHs that did not successfully merge into a SMBH 
0111, which we call IMBHs. However, despite 
the theoretical interests, it is difficult to obtain conclu- 
sive evidence for their existence. One viable strategy is 
to search for gravitational-waves produced in seed-BH 
mergers, which may become possible with the launch of 
the Laser Interferometer Space Antenna (LISA) |48ll5C| . 



B. IMBH Formation Scenarios 

In this study, we follow a previous study by Bertone 
et al. H3 and focus on two seed-BH formation scenarios 
covering the wide range of possible IMBH masses. In the 
first scenario, which we refer to as scenario A, the seed- 
BHs are remnants of the collapse of Population-Ill (or 
first generation) stars j32 ■ These stars are generally mas- 
sive due to suppression of mass-loss processes, a result of 
their very low-metal composition. As such, they are also 
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called zero-metalicity or very massive sta rs ( VMS). The 
fate of VMS have been studied by e.g. |5l],|53]: stars with 
100 < M/Mq < 250 encounter the electron-positron pair 
instability and explode in a giant nuclear-powered explo- 
sion, leaving no compact remnant, while heavier stars 
collapse completely to a BH containing at least half of 
the stellar mass |52l |. 

What we need to know is the mass-function of these 
BHs. Unfortunately, the mass-function of Population-Ill 
stars is not well known, although recent studies indicate a 
double peaked function that extends up to a few 10 3 Mm 
[53| . We base our scenario on the one proposed in [33] 
and further studied by E3, EH Hif. Interestingly, 
if BHs with masses > 10 2 M Q form at high redshifts in 
minihalos representing ~ 3<r peaks of the smoothed den- 
sity held, the resulting baryon mass fraction is found to 
be comparable to those observed for SMBHs in nearby 
galaxies As in [2l]], we conservatively consider BHs 
of mass 10 2 Mq forming in minihalos of masses larger 
than M VtCr it ~ 1O 7 M at formation redshift Zf ~ 18. 

In the second scenario, which we refer to as scenario 
B, BHs form from low angular momentum gas in proto- 
galactic disks at high redshifts, producing a population 
of seed-BHs with masses ~ 10 5 M©. We use the scenario 
proposed in |33|. which we briefly summarise. During 
the collapse of the first halos, gas cools and a pressure 
supported disk forms if the halo is massive enough to 
contain a relatively large fraction of molecular hydrogen 
(molecular hydrogen is the main coolant, see H^)- Local 
gravitational instabilities in the disk manifest themselves 
as an effective viscosity that transfers angular momen- 
tum outwards and cause an inflow of gas, in particular 
the low angular momentum tail. In halos that are both 
massive enough to contain enough hydrogen for cooling 
and do not experience mergers with other halos, the pro- 
togalactic disk can evolve uninterrupted until ultimately 
being terminated by the heating and disruption caused 
by supernovae of population-Ill stars. During this time, 
mass transfer of order ~ 1O 5 M occurs. The central mass 
may be briefly pressure supported but will ultimately col- 
lapse to a BH due to post-Newtonian instabilities. Since 
population-Ill stars have typical lifetimes of 1 — 10 Myrs, 
the halo must be avoid of mergers over many dynamical 
times. This sets a stringent lower limit on the required 
halo mass. At Zf = 12, this limit is j\I v ., r u = 10 8 M Q 
(see [33| for an expression). The masses of the BHs that 
form have a near log-normal distribution with peak mass 
= 2.3 x 10 5 Mq independent of Zf, with spread 
<Jbh = 0-9 0. 

Scenario B naturally contains many parameters other 
than Zf (such as fraction of gas cooled, lifetime of 
population-Ill stars, etc). However, uncertainties in 
these parameters largely affect Mf,h,o, and as such are 
ultimately masked by uncertainties in Zf, which affects 
the total number of sccd-BHs formed. We therefore treat 
scenario B formation to be described by one parameter, 
Zf. Now, the epoch of cosmological reionisation places 
a lower limit on zj. This is because the heating of the 



intergalactic medium and the subsequent ionisation of 
molecular hydrogen (i.e. reionisation) terminates further 
baryon cooling. Without molecular hydrogen, scenario 
B formation cannot proceed, even in the heaviest ha- 
los which satisfy M v > M v>cr it. Therefore, for scenario 
B, Zf = z re , where Zf(3 IS the redshift of reionisation 
(we assume for simplicity that seed-BH formation stops 
abruptly at reionisation). 



C. IMBH Number Density: Results from 
Numerical Studies 

Now that we have discussed IMBH formation, we will 
summarise how their number density evolves with time. 
Since in this work we are only interested in IMBHs with 
minispikes, from here on we will use 'IMBH' to im- 
ply 'an IMBH with a maintained minispike'. First, al- 
though seed-BH formation can lead to minispike forma- 
tion, various processes destroy it by varying degrees (see 
Sec. 1111 C|l . and the strength and number of minispike 
decreases with time. BH mergers, i.e. mergering with 
another IMBH or a SMBH, are the most destructive and 
we must therefore consider its effects. Before doing so 
however (in Sec. 1111 D|l we summarise basic procedures 
and results of previous numerical studies. 

The basic approach focuses on constructing a statisti- 
cal sample of halo formation histories (each called a 're- 
alisation') in the context of the hierarchical CDM model 
for structure formation, followed by computing the dy- 
namical evolution of halos and BHs within halos (see 
[H El El E3 for further details). The first step is 
to consider a virialised halo of mass M V fi at z — and 
construct a merger history, i.e. a list of all the smaller 
halos that merged together to form the final halo, as well 
as the redshifts at which the mergers occurred. The next 
step is to plant BHs in halos that satisfy seed-BH forma- 
tion conditions. That is, if a halo satisfies M v > M VyCr u 
and z > Zf, a seed-BH is planted at the earliest time 
M v > M ViCr it still holds. This point is labelled Zbh- The 
last step involves evolvin g th e halos and BHs forwards 
to z = 0, as described in [Ig,|5(j. During the last step, 
a BH that comes within a distance min(0.01?\,i r , lkpc), 
where r v i r is the virial radius of the halo, of another 
BH is considered a BH-BH 'potential pair', and is de- 
coupled from the simulation. Previous work focused on 
the gravitational- waves produced by the potential pairs' 
subsequent mergers |48j , but for our purposes we exclude 
all 'potential-pairs' from our EGB calculation. 

What necessary results can we obtain from these stud- 
ies? First, the distribution of BH formation redshifts 
Zbh is found to be exponential, i.e. using the Milky-way 
galaxy with M V £ — 10 12 ' 1 h~ 1 M ( ? ) and 200 realisations, 
|48| finds that the number of BHs formed peaks at Zf and 
decreases exponentially for higher Zbh (see their Fig. 2). 
Second, by following the dynamical evolution of BHs, 
[2l| finds that N bh ,A = 1027 ± 84 scenario A BHs and 
Nbh,B = 101 ± 22 scenario B BHs remain unmerged in a 
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Milk-way sized halo. Errors denote the la halo-to-halo 
scatter (200 realisations were performed). Finally, the 
number of 'potential pairs' formed is highest at zt when 
the BH number density is largest, and decreases as a 
power-law of (1 + z) with decreasing z (see their Fig. 3). 



III. FORMULATIONS 



dE 



(E (l + z),M bh ,z), 
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where c is the speed of light, Hq is the Hubble parameter 
now, and h(z) is the function h(z) = • v /[f2jv/(l + .z) 3 + ^A]- 
We use Eq. J2J to compute the final EGB flux. In 
the following subsections, we will deal with the physical 
quantities in the expression in detail. 



A. Extragalactic Gamma-Ray Background 



B. Dark Matter Annihilation Gamma-rays 



To calculate contributions to the EGB flux from unre- 
solved cosmological DM sources, we adopt the method- 
ology in Ullio et al. .25), but extend it for our purposes 
of including IMBH minispike enhancements. 

Let dN 1 /dE{E 1 Mbh, z) be the differential energy spec- 
trum for the number of gamma-ray photons emitted per 
unit time from a single IMBH of mass Mbh at redshift 
z, and dn/dMbh(Mbh, z) the comoving number density 
of IMBHs of mass Mbh at redshift z. From these we can 
determine the number of photons emitted in a proper vol- 
ume dV say, at redshift z in time interval dt and energy 
range (E,E + dE). Then, assuming isotropic emission, 
the corresponding number of photons d_/V 7 collected by a 
detector on Earth with an effective area dA in the red- 
shifted energy range (Eq, Eq + dEg) over a time dto, is 

d.V- = / dMb h {l + zf-^-{Mbh>z)— -3-{E,M bh ,z) 

dM h h dE 



x e 
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where we have used the fact that dEdt = (1 + z)di?o(l + 
z)~ 1 dto = dEodtg. The exponential term is an attenua- 
tion factor which accounts for the absorption of gamma- 
rays during propagation to Earth, the factor (1 + z) 3 
converts from comoving to proper IMBH number den- 
sity, and Ro comes from the metric of our cosmology as 
defined in section We define dV by radial increment 
dr and angular increment dfi as 



dV = 



(R r) 2 R Q 
{1 + z) 3 



drdfi. 



The main absorption of gamma-rays of a few GeV 
is via pair production on the extragalactic background 
light emitted by galaxies in the optical and infrared 
bands. Although at these gamma-ray energies attenua- 
tion is almost negligible, we use the form e~ z / Zma:c , where 
we approximate z max ~ 3.3(i?o/10 GeV) -08 , which is 
a parametrisation that reproduces results of |57j with 
enough accuracy. Substituting dV into Eq. QJ and 
changing the integration along the line of sight dr to 
along dz, we get for the flux 
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In order for DM particles to satisfy cosmological con- 
straints, they are expected to have a small but non-zero 
annihilation cross-section into Standard-Model particles. 
This ensures they are in chemical equilibrium in the early 
universe. By constraining their relic density by cosmo- 
logical observations, one can obtain limits on their anni- 
hilation cross-section. For DM that is a thermal relic, 
the required cross-section is 3 x 10 -26 cm 3 s -1 . How- 
ever this should be taken as an upper limit, because 
processes such as coannihilation can allow smaller cross- 
sections (see 0, Q for details) . The DM particle mass 
is constrained from below by collider experiments and 
above by theory, giving a commonly accepted range of 
m x ~ 50 GeV- lOTeV. 

Now, the flux of gamma-rays from a single IMBH min- 
ispike is quantified by the term dJ\fy/dE(E, Mbh, z) m 
Eq. J2J ■ F° r gamma-rays of DM origin, we can rewrite 
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where the factor 1/2 appears due to the fact this is an 
annihilation of identical particles j2^|, p(r) is the DM 
density profile around an IMBH ranging from radii ru m 
to r sp (we discuss these in Sec. lIII av is the annihila- 
tion cross-section times relative velocity, and dV 7 (E) / dE 
is the differential gamma-ray yield per annihilation. The 
latter can be divided into continuum and monochromatic 
(line) emissions, and can be written as [25| 



dN 7 (E) 
dE 



^2b lY n lY S(E - m x {\ - M Y /Am x )) 

Y 

+ Y^b F ^%±{E), (4) 



dE 



where the second term refers to the continuum. The 
continuum photons are produced by annihilation into the 
full set of tree-level final states F including fermions and 
gauge or Higgs bosons which generate photons on decay. 
The bulk of the photons however are produced in the 
hadronisation and decay of neutral pions (decay mode 
7T° -> 27), with a branching ratio of b F = 98.8%. The 
monochromatic emission on the other hand is the result 
of prompt annihilation into two-body states including a 
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photon, a process that is forbidden at tree-level and only 
allowed in higher order perturbation theory. Although 
subdominant, these gamma-rays have the advantage of 
producing a 'smoking gun' signature of DM annihilation, 
i.e. photons of energy E — m x for the 27 final state and 
E = m x (l — My/4rn x ) for the photon and particle Y 
final state. The parameter 6 7 y is the branching ratio 
into the respective channels, and n 7 y is the number of 
photons emitted per annihilation. 

In the current work we consider both the contin- 
uum and monochromatic gamma-ray emissions. For 
the continuum we consider the ir° branch, and use a 
conveniently parametrised form of the rest frame en- 
ergy distribution per annihilation, dN cont (E) / dE = 
(0.42/m x )e _8! 7(x 1,s + 0.00014) where x = E/m x . For 
the monochromatic emission we consider the \X ~ * ^7 
process, which exists in many supersymmetric models. 
We use for the branching ratio 627 = 10~ 3 . 



C. Dark Matter Enhancement around IMBH 

The emergence of a deeper gravitational potential due 
to the formation of a BH inevitably alters the DM halo in 
which it is formed. It has been shown that the adiabatic 
formation of a SMBH results in an enhancement of the 
nearby DM density, called a 'spike' ^4[. However, it has 
also been shown that spike formation depends on initial 
conditions |l5j |. and even formed, it is depleted by pro- 
cesses such as BH-BH mergers [lCI . dynamical processes 
such as gravitational scattering off stars 0], an d DM 
annihilation itself [TtI ]. Fortunately, minispikes around 
IMBHs may not be greatly affected by these problems. 
Firstly, we have selected IMBHs that survive without ex- 
periencing any major mergers. Secondly, these unmerged 
IMBHs are not necessarily localised in the galactic cen- 
tre, and null observations of ULXs in our galaxy suggest 
they reside in satellite halos with no significant stellar 
component. These imply that dynamical processes are 
not very significant. Thirdly, a central BH formation is a 
built-in property of scenario B, predicting a strong min- 
ispike. However, we must make clear that the precise 
likelihood of minispike formation and survival are still 
uncertain. With these in mind, we treat scenario B as 
our 'optimistic' case and scenario A our 'conservative'. 

The BH in scenario B forms in the centre of its host 
halo [33, and therefore predicts a strong minispike. To 
compute its shape we need to specify the initial DM pro- 
file. If we write the initial inner DM profile as p(r) oc r~ 7 , 
the minispike has the form [2lJ 



Ps P (r) = p(r sp ) 



-Isp 



(•5) 



where r sp is the radius of the minispike and j sp is the 
gradient, given as 



9-2 7 



(0) 



TABLE I: Representative values of halo radii for scenario B, 
for two formation redshifts. Shown are the halo virial radius 
r V i r , the scale radius r s which defines the shape of the NFW 
profile (see H3|), the 'influence' radius of the BH as de- 
fined in [oSf. and the inner radius of the minispike ru m . DM 
parameters are m x = 100 GeV and av — 3 x 10~ 26 cm 3 s -1 . 





Radii, Zf = 10.9 [pc] 


Radii, Zf = 12 [pc] 


>vir 


1300 


1100 


r s 


510 


480 


Th 


34 


33 


I sp 


6.8 


6.7 


ru m (z = 0) 


5 x 10 -3 


4 x 10" 3 



Values for 7 have been proposed by analytic fits of TV- 
body simulations, e.g. 7 = 1 by Navarro, Frenk and 
White (NFW) 58], and 7 = 1.5 by Moore et al. [59|. 
Recent simulations have produced shallower profiles (see 
[6fj| and references therein) and as well as somewhere be- 
tween 1-1.5 [^lUli^. In this work we assume a middle 
value of 7 = 1 , but keep in mind that a smaller value will 
lead to less gamma-rays. For reference, our values for 
halo radii are summarised in Table [i] 

For scenario A, while some works show that they are 
formed in the centres of their host halos I64| . others show 
fragmentation may lead to off-centre BHs 65] . Although 
motivations for a central BH are strong, to remain con- 
servative we consider an off-centre BH formation. It has 
been shown that BH formation in an uniform DM back- 
ground will form a mild r~ 3 / 2 minispike |2lLl54lp ?i|67l. 
Specifically, we consider the minispike studied in |2lLl54| . 
with p(r) = Ph(r/rh)^ 3 ^ 2 where ry L — 0.045 pc. 

In both scenario A and B minispikes, the very DM 
annihilation we are studying sets an lower limit on the 
minispike radius. Assuming DM annihilation is the 
main process by which the inner density decreases, the 
DM number density n x obeys the evolution equation 



h x (r,t) 



-avn x (r,t). Solving this yields the solution 



n x (r,t) 



c(r,t f ) 



l + n x (r, t f )av{t~t f Y 



(7) 



where t — tf is the time elapsed since BH formation. The 
upper limit is of order m x /av(t — tf), and we define ru m 
as the radius at which the following holds 



Psp (flim ) 



crv(t-tf) 



= Pli 



(8) 



where m x is the DM mass. For common DM parameters, 
has grown to ~ 10~ 3 pc by z = 0. 



n 



D. Modelling the IMBH Number Density 

Now we parametrise the decreasing IMBH number den- 
sity, using results of previous numerical studies. To deter- 
mine the initial (i.e. at Zf) IMBH number density, we con- 
sider a delta- function seed-BH formation occurring at Zf, 
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and plant BHs in halos that satisfy M v > M v ^ cr i t . This 
simplified picture neglects BHs that would have formed 
earlier. That is, for some halos we plant one BH where 
in fact there would be two (or more). However, as |4S[ 
shows, the number of BHs formed in z > Zf decreases 
exponentially. This, coupled to the fact that massive 
halos become increasingly unlikely for higher redshifts, 
the effect of ignoring them is minimal. We also stress 
that this picture in no way overestimates EGB contribu- 
tions, because (i) we have underestimated the initial BH 
density, and (ii) we have neglected all DM annihilations 
before Zf. Moreover, any underestimation is masked by 
uncertainties caused by Zf. 
The required calculation is 



n bh (z f ) 



d.72 



where n b h is the comoving number density of IMBHs that 
have not experienced any mergers, and dn/dM v is the 
halo mass-function. For the mass-function we use one 
postulated by Press-Schechter theory j6(|, 



dn 
dM 



_ Po f/ ' 



dlog^ 



'dlogM' (10) 

where p is the comoving matter background density, 
p = pcf^M, with p c the critical density. The parame- 
ter v = 5 sc (z)/a(M) is defined as the ratio between the 
critical linear fractional overdensity required for vitalisa- 
tion over a(M) the present rms linear density fluctuation 
in spheres containing a mean mass M. For the multiplic- 
ity function f(v) we use the ellipsoidal collapse model 
[7C| , normalised to results of A-body simulations of the 
Virgo Consortium [7l]], as was done in plj . 

Now, we assume that the number of IMBHs present 
and the number of 'potential pairs' formed are propor- 
tional, and fit a power-law redshift dependency to n b h, 



n bh (z) = n bh (z f ) 



l + z 



1 



z f 



(11) 



where f3 is a free parameter. We can derive representative 
values for (3 by fitting n b h (0) to numerical results by [2lJ . 
From their value of N b h t A {Nbh,B)-, we interpolate the 
number of unmerged IMBHs in other galaxies by their 
halo masses, and determine n^O) as 



M„, ( 



M v 



n bh (0) = N bh dM v in12-1 " -— (0), 



10 12 - 1 h~ 1 M G 







dn 



dM, 



(12) 

where a is some constant which we assume to be 1. This 
gives (3 ~ 0.2 (0.3) for scenarios A (scenario B). Although 
in reality we do not expect a to be exactly 1 (e.g. we ex- 
pect less unmerged IMBHs in older elliptical galaxies), 
we find that for reasonable deviations from a — 1 our 
estimates of (3 do not appreciably change. Moreover, un- 
certainties caused by deviations in a are far smaller than 
those caused by Zf. 

Lastly, we conclude this section with the parameters 
of scenarios A and B in Table HH 



TABLE II: Scenario A and B parameters, including formation 
redshift Zf, the minimum halo mass required for formation to 
occur M VtCr it, the number density of IMBH n b h at z = Zf 
determined from Eq. @ , the number of IMBHs residing in a 
Milky- Way size galactic halo Nbh [2l|. n bb/ at z = derived 
from Eq. 1121 . and /3 derived from our fitting. 



z f 

M bh [M Q ] 

M v , cri t(Zf) [Mq] 

n bh (z f ) [Mpc" 3 ] 
n bh (0) [Mpc" 3 ] 



a z re taken from WMAP 3rd year results |3 

IV. CALCULATION AND RESULTS 

A. Contribution to the EGB 

Substituting Eq. J2J into Eq. © we find the gamma- 
ray flux observed at Earth as 



Scenario A 


Scenario B 


18 


10.9+^ ° 


10 2 


10 5 ' 


4 x 10 B 


10 8 


23 


2.5 


1027 ± 84 


101 ± 22 


12 


1.1 


0.2 


0.3 



d$ 7 av c Po 



dz A2iz)dN : {E,) c-^ 



dE 8vr H ra\ J " h{z) dE 
where E 1 = Eq(1 + z) and A 2 (z) is defined 
dn 



A 2 (z) = i 
Po 



dM, 



hi, 



dM, 



hi, 



(13) 



(M bh ,z) / p\r)d 6 r. (14) 



Note the disappearance of the (1 + z) 3 term in Eq. i|13|) 
compared to [25|]. The reason for this is that while the 
DM density in the halo is a function of z, the density in 
the minispike around an IMBH is only a function of Zf] 
the z dependence is taken into account instead by ru m . 

For scenario B, since the BH mass distribution is near 
log-normal, and the enhancement due to the minispike 
scales linearly with the BH mass, we find that with 
very good accuracy we can approximate all the BH mass 
to have the peak value M^o- Hence, we substitute 
dn/dM bh (M bh ,z) = n bh (z)S(M bh - M bh>0 ), which gives 
the final form we will use, 



\ z ) = ZJ2 n bh{z) 

Po 



i 2 d 3 7 



(15) 



In Fig. n we show the redshift dependence of the en- 
hancement function A(z) 2 / 'h(z) for scenario B, for three 
values of z/ = z re = 8, 10, and 12. The strongest en- 
hancement occurs at Zf, partly because of the highest 
IMBH number density, but mainly because of the pres- 
ence of the sharpest minispike. This displays the merits 
of considering the EGB, because high redshift enhance- 
ments are most easily observed as a diffuse background. 
On the other hand, we must therefore carefully consider 
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FIG. 1: Enhancement factor A 2 as a function of redshift, for 
scenario B IMBH minispikes. Three values of Zf = z re are 
plotted, 8, 10, and 12. The peak enhancement at zj is largely 
due to the fact that the minispike is sharpest just after it is 
formed. This demonstrates the effectiveness of the minispike 
in enhancing contributions to the EGB. 

the IMBH number density at Zf. However, a major un- 
certainty in the hierarchical formation scenario of SMBHs 
is that the 'occupation number', that is, the fraction of 
galaxies containing a seed-BH at high redshift, is highly 
model dependent. In other words, contributions to the 
EGB flux can potentially constrain seed-BH models. To 
this effect, our calculation for the initial density is a con- 
servative value (see Sec. IIIID|) . 

In Fig. |21 we show contributions to the EGB from 
DM annihilation into continuous gamma-rays (we discuss 
monochromatic emission later). We use DM parameters 
m x = 100 GeV and av = 3 x 10~ 26 cm 3 s" 1 . The values 
of Zf and (3 in Table llll were used. We do not expect low- 
redshift IMBHs to be resolved; it has been shown that 
GLAST will resolve minispikes in our galaxy and possi- 
bly Andromeda, but not further [2l]]. Moreover, we find 
that sources within z — 0.01 contribute less than a tenth 
of the total EGB contribution. We therefore show the 
flux, integrated safely from z = to Zf. We find that 
minispikes increase the gamma-ray flux from DM halos 
by 1 — 3 orders. In particular, scenario B can give fluxes 
that are of the order of current observations. Scenario A 
fluxes are two orders smaller, but we stress that scenario 
A is a conservative case, using the smallest IMBH mass 
and a mild r~ 3 / 2 minispike. The 'host halo only' is a pre- 
diction without any spikes nor minispikes, and therefore 
acts as a minimal prediction. The dashed curves indicate 
uncertainties caused by la scatter in: Nbh,A for scenario 
A, and N^^b and z re combined for scenario B. 

In Fig. |3| we show how EGB contribution from scenario 



FIG. 2: Contributions to the EGB flux, from scenario A and 
scenario B IMBH minispikes. Also shown are EGRET data 
and predictions of our minimal host halo only scenario (i.e. no 
spikes and minispikes). We see that minispikes increase EGB 
contributions by 1 - 3 orders in magnitude. The la scatter 
in Nbh,A |2ll| is shown for scenario A. For scenario B, the la 
scatter in Nbh b as well ELS Zre £1X6 shown combined. In all 
calculations, zj and /3 shown in Table ITT1 are used, with DM 
parameters m x = 100 GeV and av — 3 x 10~ 26 cm 3 s _1 . Note 
that scenario A is our conservative case; a BH of mass 1O 2 M0 
with a mild r~ 3 ^ 2 minispike. 



B depends on the free parameter (3. As before, we use 
z re = 10.9, m x = 100 GeV, and av = 3 x 10~ 26 cm 3 s _1 . 
We show the la scatter in N b h y B using vertical dashed 
lines. We find that our fitting is within la error of 
EGRET observations. 



B. Constraining DM Parameters 

Instead of performing a complete scan over DM param- 
eters space, we follow Bertone et al. [2lJ and consider two 
discrete cases. In addition to our previously chosen set 
m x — 100 GeV and av = 3 x 10 _26 cm 3 s _1 , we define 
a new set m x — 1 TeV and av = 10 _29 cm 3 s _1 . Using 
scenario B with z re — 10.9 and (3 — 0.3, we plot in Fig.Q] 
the predicted EGB flux for both sets. 

Now in Fig. [3] we show an advantage, from the point 
of view of constraining DM parameters, of considering 
IMBH minispike enhancements. We show contributions 
to the EGB at an energy of 1 GeV, against the DM anni- 
hilation cross-section. As before, we use z re — 10.9 and 
j3 = 0.3, and assume for now m x — 100 GeV. The ad- 
vantage of minispikes is observed in the gradient of the 
plotted line. Although one would naively expect gamma- 
ray fluxes to scale as av /m 2 , the presence of a minispike 
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FIG. 3: Plot showing how EGB contributions from scenario 
B IMBH minispikes varies due /3, the parameter of our IMBH 
number density fitting. Plotted is the flux at an energy of 
1 GeV, where the contribution is greatest. The range of j3 
determined from the la scatter in Nbh,B is shown by the ver- 
tical dashed lines. EGRET observations are shown by the 
horizontal rectangle. 



FIG. 4: Contributions to the EGB from scenario B IMBH 
minispikes, for two sets of DM parameters, the first with m x = 
100 GeV and av = 3 x 10~ 26 cm 3 s~ x and the second with 
m x = ITeV and av — 10 _29 cm 3 s _1 . Errors include those 
from Nbh,B and z re . 



alters this dependence. This is due to the fact that the 
dominant term in the minispike enhancement [the inte- 
gral in Eq. (|15[) ] is given by the expression ~ Pu m r iim' 
which brings the gamma-ray flux scaling as (av) 2 ^ 7 m^/ 7 , 
for minispikes growing out of a 7 = 1 profile. Physi- 
cally speaking, a smaller cross-section works to maintain 
a denser minispike, which compensates for the decrease 
in flux due to the smaller cross-section. One can say this 
weak dependence on DM parameters make minispikes 
particularly suited for DM detection. 

To determine the potential of the weak dependence 
on DM parameters, we consider gamma-ray observations 
by GLAST, which is expected to have more than an or- 
der better point source sensitivity than EGRET. With 
its launch, many gamma-ray sources that could not have 
been resolved until now will be detected, and taking these 
into account, a smaller EGB flux is expected. The most 
widely considered candidate for the dominant EGB con- 
tributor is unresolved blazars^.e. a beamed population 
of active galactic nuclei (see [73 and references therein) , 
and the fraction of these blazars that can be removed 
with GLAST depends on their luminosity function. The 
latest calculation predicts that the resolvable fraction is 
around 20% If the EGB flux is entirely due to 

blazars, then the EGB will be reduced by the same frac- 
tion, and hence the sensitivity to av will be improve d by 
a factor of 2. We should, however, keep in mind that [73] 
concluded that blazars cannot fully explain the EGB flux, 
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FIG. 5: Gamma-ray flux at E = 1 GeV from scenario B 
IMBHs, plotted against the annihilation cross-section av. 
The dashed lines indicate combined error in Nf,h,B and z re - 
EGB observations by EGRET and av of 'natural' thermal 
neutralinos are shown. The weak dependency on av works 
positively for DM detection. 
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but only 25-50%. Therefore, the remaining 50-75% may 
be due to other astrophysical objects of either known or 
unknown origin, which may or may not be resolved by 
GLAST. This indicates that a significant fraction may 
still be resolved with GLAST, depending on the prop- 
erty of this additional contributor. In an optimistic case, 
where a fraction of 0.75 + 0.25 x 0.2 = 0.8 can be resolved 
with GLAST, the cross-section sensitivity will be around 
3 x 10~ 29 cm~ 2 s -1 , depending on values of Zf and N b h 
(see error bars in Fig. [SJ. The sensitivity will also de- 
pend on M bh ; if halve M bh , the sensitivity decreases to 
~ 10~ 28 cm~ 2 s _1 . These sensitivities are so small that 
no other experiments can compete in the next decades. 



C. Line Gamma- Rays 

In the previous subsection we discussed constraining 
DM parameters using contributions to the EGB from DM 
annihilation into continuous photons, and highlighted its 
potential using the future GLAST mission. However, 
the continuum component lacks a distinguishing signa- 
ture to separate it from other sources, such as unresolved 
blazars. Although the spectrum of other sources de- 
creases rapidly at high energies, this is nonetheless a dif- 
ficulty. Ideally, we would like to identify DM annihilation 
without these complications. The monochromatic com- 
ponent of DM annihilation provides a means of achieving 
this. As discussed earlier, the monochromatic emission 
provides a 'smoking gun' signature of DM annihilation 
due to its energy at the DM mass. Here, we present 
the monochromatic spectrum from the XX ~ ¥ ^1 pro- 
cess with a fixed branching ratio of 627 = 10~ 3 . Our 
result is shown in Fig. for two sets of DM parameters, 
m x — 100 GeV and av = 10~ 26 cm 3 s _1 (note this is not 
excluded by the continuum component), and m x = 1 TeV 
and av = 10 _29 cm 3 s~ 1 . Again, we consider scenario B 
IMBH minispikes with z re — 10.9 and j3 = 0.3. 

The spectral shape of Fig.Elarises because of distortion 
due to cosmological redshift and absorption of gamma- 
rays during propagation. The larger peak at Eq — m x 
is characteristic of DM annihilation, and if detected pro- 
vides very convincing evidence for DM. The smaller peak 
is characteristic of IMBH minispikes; it is due to the 
strong enhancement factor at Zf, as shown in Fig. ^ If 
detected, this second peak identifies the presence of a 
high enhancement factor at high redshift, supporting the 
case of seed-BHs and minispikes. 

In the past, line gamma-rays had not been intensively 
considered, because the EGB spectrum taken by EGRET 
did not reach up to the required O(100) GeV energies. 
With GLAST however, the energy window extends up to 
300 GeV (and with better energy resolution), and line 
gamma-rays could be a serious candidate of the first DM 
detection. Furthermore, as shown in Fig. the min- 
ispike around IMBHs provide such high EGB line fluxes 
that it might give better evidence than the continuum 
flux. The characteristic spectral feature, combined with 




10" a 0.1 1 10 100 1000 
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FIG. 6: Spectral signature in the EGB due to DM annihilation 
into monochromatic photons in scenario B IMBH minispikes. 
The 27 branch is considered here, with a fixed branching ratio 
of 10~ 3 . Two DM parameters sets are shown, m x — 100 GeV 
and av = 10 -26 cm 3 s _1 (this is not excluded by the contin- 
uum component), and m x = ITeV and av — 10~ 29 cm 3 s _1 . 
We used z re = 10.9 and /3 = 0.3. 

good energy resolution of GLAST, work quite positively 
for the detection of line gamma-rays, even if gamma-rays 
from other astrophysical sources give considerable con- 
tribution at the same energy. 

V. DISCUSSION AND CONCLUSIONS 

We have studied contributions to the EGB from DM 
annihilation in minispikes around IMBHs. Our results 
are plotted in Fig. fusing DM parameters m x = 100 GeV 
and av = 3 x 10~ 26 cm 3 s _1 . We find that a consideration 
of minispikes increase the contribution by 1 — 3 orders, so 
that in optimistic scenarios the predicted gamma-ray flux 
may reach current EGB values. The EGB can therefore 
potentially be used to constrain DM parameters and/or 
IMBH scenarios, particularly when better EGB observa- 
tions are taken by GLAST. In our work, we considered 
two IMBH formation scenarios, scenario A being rem- 
nants of Population-Ill stars [32] , and scenario B being 
formed in the centres of protogalactic disks 33]. Sce- 
nario A is our conservative case (mass 10 2 M Q with a 
mild r -3 / 2 minispike) while scenario B is our optimistic 
(mass 10 5 Mq with a r~ 7 / 3 minispike). 

We also showed that DM annihilation into monochro- 
matic gamma-rays may be a serious contender for indi- 
rect DM detection. Our result, for scenario B IMBHs 
and two sets of DM parameters, is shown in Fig. Note 
that these parameters were chosen so that it is not ex- 
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eluded by the continuum component. The higher energy 
peak, at an energy equivalent to the DM mass, is within 
GLAST's potential observation range, and can provide a 
potential 'smoking-gun' signature in the EGB. 

Compared to SMBH spikes, IMBH minispikes have the 
disadvantage that it is smaller, and it grows out of less 
dense DM profiles. Also, the survival probability of spikes 
and minispikes are still uncertain. However, there are 
still advantages to using IMBH minispikes. First, it has 
been shown that observations of our galaxy's centre con- 
strain the strength of spikes to such a degree that DM 
annihilation in spikes cannot significantly contribute to 
the EGB [2^|. Second, although survival probabilities 
are uncertain, we can select IMBHs that are likely to 
have maintained their minispikes. Our selection involves 
choosing IMBHs that have not experienced any mergers, 
on the grounds that mergers strongly destroy DM en- 
hancements. Additionally, such unmerged IMBHs are 
likely to reside in the outskirts of galactic halos, where it 
is not affected by dynamical process that are also known 
to deplete DM enhancements. 

It must also be added that minispike formation re- 
quires some initial conditions to be met, including an 
adiabatic and symmetric BH formation, formation of a 
BH in the centre of its host halo, and very cold initial DM 
orbits near the halo centre. These conditions are gener- 
ally supported by the collisionlcss nature of particle DM, 
and adiabaticity is satisfied when one compares the BH 
formation time-scale to the dynamical time-scale at some 
relevant distance |2l|. However, the collapse and accre- 
tion processes during BH formation are not well known, 
and are likely to be complex and far from symmetric. 
Also, we have neglected the effects of seed-BHs born with 
enough kick-recoil velocities to be expelled out of their 
host halos. Although a detailed study is beyond the scope 
of this paper, minispike formation for such seeds may be 
surpressed. 



An advantage of DM detection using minispikes is the 
fact that the gamma-ray flux is weakly dependent on 
DM parameters. This is because the smaller the av, the 
longer the minispike remains sharp. This fact compen- 
sates for the lack in flux, bringing the flux proportional 
to (ff!)) 2 / 7 mf. As a consequence, GLAST may be able 
to probe down to av ~ 10 _29 cm 3 s _1 if we optimisti- 
cally assume scenario B IMBH minispikes, and that 80% 
of the EGB is resolvable by GLAST. This is far smaller 
than any other experiment, and excludes the allowed av 
region due to the standard thermal relic SUSY DM sce- 
nario. However, evidence concerning IMBH properties 
and EGB composition are required before setting such 
constraints. In this respect, it has recently been shown 
that the EGB power spectrum can be used to discrimi- 
nate DM contributions from other sources [3l| . 

All our calculations can be applied to other numer- 
ous DM candidates by simply substituting the appro- 
priate differential gamma-ray spectrum per annihilation 
Eq. J3J|. A detailed search over particle DM candidates 
using EGB contributions should become possible with 
the launch of GLAST, and these should then be cross- 
correlated with constraints from other potential DM sig- 
nals. Finally, we have used EGB contributions to reveal 
DM properties, but we stress that the fact that scenario 
B predicts gamma-ray fluxes that are two orders greater 
than scenario A may also shed light on IMBH and SMBH 
models. 
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